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Abstract 

We show that any solution of the SU (2) Skyrme model can be used to give a topologically 
trivial solution of the SU(4) one. In addition, we extend the method introduced in [|l] 
and use harmonic maps from 5 12 to CP^" 1 -* to construct low energy configurations of the 
SU (N) Skyrme models. We show that one of such maps gives an exact, topologically trivial, 
solution of the SU(3) model. We study various properties of these maps and show that, 
in general, their energies are only marginally higher than the energies of the corresponding 
SU(2) embeddings. Moreover, we show that the baryon (and energy) densities of the SU(3) 
configurations with baryon number B = 2 — 4 are more symmetrical than their SU (2) 
analogues. We also present the baryon densities for the B — 5 and B — 6 configurations and 
discuss their symmetries. 



1 INTRODUCTION 



The Skyrme model is now well established as an effective classical theory used to describe 
nuclei 0, |3j for which the field, which describes pions, is valued in SU(N). To have finite 
energy configurations, one must require that the field U(x, t) goes to a constant matrix, 
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say J, at spatial infinity: U — > I as \x\ — > oo. This effectively compactifies the three 
dimensional Euclidean space into S 3 and hence implies that the field configurations of the 
Skyrme model can be considered as mappings from S 3 into SU(N). 

In terms of the algebra valued currents d^UU -1 , the Lagrangian is, 



where U(x,t) is an SU(N) valued scalar field, F ~ 189 MeV is the pion decay constant 
and a is a dimensional constant. 

The last term, called the Skyrme term, stabilises the solitons and, in addition, in- 
troduces small interaction forces between them. Their nature depends on the relative 
orientation of skyrmions in their internal space. If we want to study interactions of physi- 
cal mesons we have to introduce further terms which are responsible for the meson masses, 
ie terms of the form 

C m = ^ Ml tr (U- 1 + U-2I). (2) 

Such terms play a more significant role in lower dimensions since in (2+1) dimensions their 
presence, together with the Skyrme term, is required to stabilise the solitons. However, 
in (3+1) dimensions the solitons are stable even when M v = 0. In what follows, we will 
set M n = except in section 5. 

In this paper, we concentrate our attention on studying the static properties of the 
model and we consider U(x) fields, which are stationary points of the energy functional. 
After rescaling the space-time coordinates M, E| x — > 2x/aF and defining = M T 2/aF, 
we can write the static energy corresponding to the lagrangian (||) and @ as 

E = — [ {--^(diUU- 1 ) 2 --tv 
AaJB?\ 2 V • J 16 

(3) 

From now on, we will take F/4a = 1/127T 2 so that the energy is expressed in the same 
units as the baryon number. In this case, the static fields U obey the equation 

d t (diUU- 1 - l - [djUU- 1 , [djUU^tdiUU- 1 ]]) -^(U- U- 1 ) = 0. (4) 

As the third homotopy class of SU(N) is Z every field configuration is characterised by 
an integer: 

B= 2^L £l]k tr ( diU U ~ ld i U U ~ ld * U £S > ( 5 ) 
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= — r tr fef/ ^f/- 1 ) + -— tr d^UU'\d u UU 
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,djUU 



-i 



^tr(C7- 



+ U - 2 I) )£ 



x. 



which is to be interpreted as the baryon number [[| |J; therefore, the lowest energy state 
in the B = 1 sector can be identified with the (classical) nucleon. 

So far most of the studies involving the Skyrme model have concentrated on the SU{2) 
version of the model and its embeddings into SU(N). The simplest nontrivial classical 
solution involves a single skyrmion (B = 1) and has already been discussed by Skyrme 
0. The energy density of this solution is radially symmetric and, as a result, using the 
so-called hedgehog ansatz one can reduce (Hp to an ordinary differential equation which 
can then be solved numerically. 

Many solutions with B > 1 have also been computed numerically and in all cases the 
solutions are very symmetrical (cf. Battye et al. || and references therein). The energy 
density of the two skyrmion solution forms a torus, while the energy density of the B = 3 
solution has the symmetry of a tetrahedron. For larger B the solutions describe semi- 
radially symmetric structures in which skyrmions break up into connected parts which 
are all located on a spherical hollow shell and, as was shown in ]l| , the positions of these 
skyrmionic parts on S 2 are very symmetrical. 

Recently, Houghton et al. |]J have showed that by using rational maps from S 2 to 
S 2 one can easily construct field configurations for the SU(2) model which are close to 
being solutions of the model: they have energies slightly higher than the energies of the 
exact solutions found numerically but the symmetries of the baryon and energy densities 
are the same. When these configurations are used as initial conditions in a relaxation 
program, the fields do not change much as they evolve towards the exact solutions. 

All this work has involved the SU(2) skyrmions; however, so far, very little has been 
done for the SU(N) model when N > 2. An interesting question then arises as to 
whether there are any finite energy solutions of the SU(N) (N > 2) model which are not 
embeddings of the SU(2) model and, if they exist, whether they have lower energies than 
their SU{2) counterparts. 

The first example of such a non- embedding configuration for a higher group was the 
SO (3) soliton, which corresponds to a bound system of two skyrmions, and which was 
found using the chiral field ansatz by Balachandran et al. [fjj]. Another configuration, 
with a large SU(3) strangeness content, was found by Kopeliovich et al. |§. However, all 
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other known skyrmion configurations seem to have been the embeddings of the solutions 
of the SU{2) model. 



2 SU(2) EMBEDDINGS 



As we have said before, any solution of the SU(N') model is automatically a solution 
of any SU(N) model as long as N' < N; simply by completing the entries of the larger 
matrix with l's along the diagonal and O's off diagonal. The energy, the baryon number 
and all other properties are unchanged by this operation and so such embeddings have 
the same properties as the original fields. 

However, there exist further, less obvious embeddings in which new fields have different 
properties from the original ones. In particular, one can show that any solution of the 
SU(2) model generates a solution of the S77(4) one. 

A special feature of the SU (2) field is that it can be written as 

U = f r , (6) 

where f = (1, ia±, ia 2 , i<7 3 ) and the cr's stand for the Pauli matrices. The unitarity of U 
requires that • = 1; and the energy density of this SU{2) field is 

S 2 = (34- ,9,0) + X - [(dif ■ <9,0) 2 - (<9,0 • d^f] . (7) 

Moreover, the equations of motion which follow from (HI) are 



a + <Pu {<Pj ■ <Pj) + <Pi Wij ■ <Pj) - <Pij {<Pi ■ <Pj) - <Pj {<Pii ■ <Pj, 



-0 (0 • 0jj) (0j • 0,') - (0-0^) (0, • (f)j) 







where 0j = <9,0 and 0,j = 9^0. 



Notice now that we can construct an SU(A) field out of any S* 3 field by taking 

U=U (7-20®^), (9) 

where U is a constant matrix. Note that det U = — det U and so by choosing U to 
be a constant matrix of determinant -1 we see that U is unitary. In this case, U^ 1 = 
(I - 2 <g> f ) Uq 1 and the £/" 's cancel in (|). 



To derive the equation that the field must satisfy so that (H) is a solution of 
with = 0, we note that the condition 0-0 = 1 gives 

djUU- 1 = 2<9 i 0®0 t - 20® o^, (10) 

and so (|]) becomes 

2di [0 ® 6^ (^-0 • <9i0) -0 ® 9^ (o>,-0 • ^-0) +^0 ® f (dj-0 • 0^-0) -a,-0 ® f (^0 • ^-0)" 

+di(d i 0®0 t -0®<9 i t ) =0. 

(11) 

One can easily show that the equations ( [TT]) and @ are equivalent, implying that 
satisfies the equation of the SU(2) Skyrme model (|TTD. We can thus conclude that any 
solution of the SU(2) model can be transformed into a solution of the SU(4) model by 
the embedding (g). 

The solutions obtained this way are topologically trivial since their baryon density 
vanishes identically. Moreover, their energy density is 

£ 3 = A (84 ■ 84) + 2 [(84 ■ d4f - (84 • <9;0) 2 ] , (12) 
ie it is four times larger than the corresponding energy of the original SU(2) field (0). 

This suggests that these particular SU (4) solutions may be interpreted as states corre- 
sponding to 2B skyrmions and 2B anti-skyrmions, where B is the baryon number of the 
original SU(2) solution. Incidentally, a similar situation arises in 2-dimensions where any 
B solitonic solution of the CP 1 model gives a topologically trivial solution of the CP 2 
model which can be interpreted as a bound state of 2B solitons and 2B anti-solitons Q . 



3 HARMONIC MAPS ANSATZ 



Recently, Houghton et al. [jl| exploited the similarity of the energy densities of the multi- 
skyrmion solutions to those of the SU(2) BPS monopoles and presented a new ansatz 
for constructing SU(2) multi-skyrmion fields; based on rational maps of the two dimen- 
sional sphere S 2 . Namely, they showed that solutions of the Skyrme model can be well 
approximated by the expressions of the form 

U(r,9,(j)) = exp(ig(r) h • a), (13) 
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where (r, 9, 0) are the usual polar coordinates on R 3 , and 

n = TT ^(2K( J R), 29f(i2), 1 - \R\ 2 ), (14) 

where R are rational functions of £ = tan(#/2)e l<?!> and where g(r) is a real function satis- 
fying the boundary conditions: g(0) = ir and g(oo) = 0. In other words, the configuration 
(|TB|) involves a radial profile function g{r) and a rational map from the two dimensional 
sphere of radius r which can be identified with a sphere centered at the origin, in R 3 , into 
a S 2 submanifold of SU(2) = S 3 . Moreover, it is easy to check that the baryon number 
B is given by the degree of the rational map n. 

To determine g and h one must insert (|1^) into (|j) and minimise the energy. It turns 
out that, for this minimum, n must be a rational map with a large discrete symmetry 
and that g(r) satisfies an ordinary differential equation. 

In this section we show that this idea of Houghton et al. can be generalised to SU(N). 
Using the polar coordinates (r, 9, <fi) in R 3 , our generalisation of Houghton et al.'s ansatz 
is to consider U of the form 

U{r,9,<f>) = e ^(r)(P-W) 

= e~ 2i9{r)/N (/ + (e 2i9 - 1)P) , (15) 

where P is a N x N hermitian projector which depends only on the angular variables 
(0,4>) and g(r) is the radial profile function. Note that, the matrix P can be thought of 
as a mapping from S 2 into CP^ N ~ l \ Hence it is convenient, rather than using the polar 
coordinates, to map the sphere onto the complex plane via a stereographic projection 
and, instead of 9 and </>, use the complex coordinate £ and its conjugate. Thus, P can be 
written as 

V <g> 

p (v) = (16) 

where V is a iV component complex vector (dependent on £ and £). 

For ([TBI ) to be well-defined at the origin, like (0), the radial profile function g(r) has 
to satisfy g(0) = n while the boundary value U — > I at r = oo requires that g"(oo) = 0. 
An attractive feature of the ansatz QI5D is that it leads to a simple expression for the 
energy density which can be successively minimized with respect to the parameters of the 
projector P and then with respect to the shape of the profile function g(r). This is then 
expected to give good approximations to multi-skyrmion field configurations. 
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Moreover, we will show that this method not only allows us to find such field config- 
urations but also gives us an exact non-topological solution of the £77(3) Skyrme model. 
We will also present some upper bounds on the energy of some multi-skyrmion field con- 
figurations in the SU(N) model (with radially symmetric energy density distribution). In 
what follows, we restrict our attention to the case = 0. 

To find the exact solution of the SU(3) model we put (|T5|) into (|J) and obtain 



-3 d r (9rT 2 ) (I - P) - d r (g r \A\ 2 ) (1 + |£| 2 ) 2 ([P b PPz\ + [P^PPi 



A\A\ 2 

' 16r 4 ' 
\A\ 



16r 4 



i + iei 2 ) 2 {<% ((i + [p*, m) -o, ((i + \^m, m)} 



2\2 



where A = e 2,19 — 1. 



Moreover, the energy (^) simplifies to 
1 



E = — J dr ( A N g 2 r 2 + 2J\f sin 2 g(l + g 2 ) +1 



sin 4 g " 



where 



(17) 



18) 



An = jj{N-l), 
U = ^/^^tr(|a 5 P| 2 ) 



1 = ^/^rfe(i + iei 2 ) 2 tr([^p,%p] 



(19) 
(20) 
(21) 



As the integrals M and X in ([18|) are independent of r, we can minimise ([H]) by first 
minimising M and X as functions of P and then with respect to the profile function g. 



However, since is the expression for the energy of the two dimensional Euclidean 
CP 2 sigma model, all classical solutions contain the so-called self-dual solutions, instan- 
tons or holomorphic maps from £ 2 into CP^ N ~ X \ first given in [11], which are given by 
the projector P of the form ([II]) with V — /(£)■ I n this case, the energy M is given by 
the degree of /, ie the degree of the highest order polynomial in £ among the components 
of / after all their common factors have been canceled. 
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By a Bogomolny-type argument it can be shown that 



(2 \ " 
g r r^A^ + VT^^-j + 2iVsin 2 #(l + g r f - 2g r sin 2 g ^2N + ^A N X^ 



> - \2N + \IA N Z ) 



(22) 



Finally, the baryon number for this ansatz is given by 



B 



d£ d£ tr [P [d ( P, %P] J / dr sin 2 g g 



= ^JdZdZti(p[dtP,dzP]), 

which is the topological charge of the two-dimensional CP^ N ^ sigma model. 



(23) 



In the next two sections we will show that this ansatz gives us interesting low energy 
field configurations of the SU(N) Skyrme model which are not the 577(2) embeddings. 
To minimise E we will, first of all, fix the baryon number J\f — B of the configurations 
we are interested in. We will then minimise I over all maps of degree B and then derive 



a second order differential equation for g by minimising the energy (18) treating M and 
X as parameters. 



4 5*7(3) EXACT SOLUTION 



When the projector P is analytic, ie is of the form 

P o = P{f) = Tjl, (24) 

where / is a holomorphic vector (ie whose entries are functions only of £) then it satisfies 
the equation 

P %Po = 0, %P Po = <%P , (25) 

ie the self-dual equations of the two dimensional CP^ N ~^ sigma models ||. 

Following [[Kj, we define the operator P + by its action on any vector v G as 

P+v = d(V - !^1M, (26) 
\v\ 

and then define further vectors P+v by induction: P+v = P + (P+~ 1 v) . 



To proceed further we note the following useful properties of P^f when / is holomor- 
phic: 

(P*tfP l + f = Q, k^l, (27) 

f k \ k \P+f\ 2 ( P+ l f \ P+f 

<% \ P +f ) = ~ P + f |pfc-lj| 2 ' ( |pfe-ly|2 J = |pfc-lj|2' ( 28 ) 

These properties either follow directly from the definition of the P + operator or are very 
easy to prove M. 

It is also convenient to define projectors corresponding to the family of P*/ vectors 
as follows: 

P = P(/), P 1 = P(P + f), P = P(P;/). (29) 



Taking P = Pj, for given and using the above properties we observe that all the 
terms in (|I7|), except the first one, can be gathered into one term if and only if 

l^ff = K (30) 
l/l 2 (1 + KI 2 ) 2 ' 1 ) 

where K is a constant. Moreover, for the SU(2) case, the projectors Pq and Pi satisfy the 

relation 

Po + Pi = /, (31) 

and for / = (1,£) all the terms in ([17]) are proportional to one common matrix thus 
giving a second order differential equation for the profile function g. This means that the 
Skyrme field (|15D , in the case when g satisfies its equation, is an exact solution of the 
equation (|17|). A little thought shows that this is the well known hedgehog solution. 



Unfortunately this discussion does not generalise to higher SU(N) groups. However, 
we note that for the SU(3) model, if we take P = Pi and use the fact that P0+P1+P2 = /, 
all the matrix terms in equation ([17]) become proportional to each other leading to a second 
order differential equation for the profile function, if and only if 

\m? + s ^4^, (32) 



\p+f\ 2 1/1 2 ' (1 + ia) 2 

where K, is a constant. This last condition is satisfied if 



f = (i, ey. 
9 



(33) 



Thus, by taking P = P\ for / of the form (p3[), and requiring g to satisfy the equation 



(\ sin 2 o\ 2 q r sin2o / „ sin 2 o\ . . 

g rr - + 2— * +-^ + ^ [g 2 r -l ^ =0, (34) 



we see that ( |i~5|) is an exact solution of the SU(3) model. 

For this solution, the parameters in the energy density can be evaluated analytically; 
we find 

a n = \, AT = 4, X = 4, (35) 
and the total energy is E = 3.861. 

To understand what this solution corresponds to we calculate the topological charge 
of this configuration and find 

B =y^(!W J if)' (36) 

which due to the conditions (p7|), (f2~8|) and (|32"D is identically zero. 



Although the baryon density is identically zero the solution itself is nontrivial. This 
follows from the fact that the CP 2 sigma model harmonic map P\ corresponds to a mixture 
of two solitons and two anti-solitons. Thus it seems reasonable to interpret this solution 
as describing a bound state of two skyrmions and two anti-skyrmions and as such to be 
unstable, ie correspond to a saddle point of the energy. However, let us emphasize, once 
again, that this field configuration is a genuine solution of the SU(3) Skyrme model. 

It is easy to see that this new field configuration has an energy density distribution 
shaped like a shell (ie is radially symmetric). To see this note that for this solution, 
tr(|<%P| 2 ) and tr ([<9 5 P, <%P] 2 ) which appear in (^) and (j2TJ), are proportional to (1 + 



|£| 2 ) 2 and (1 + |£| 2 ) 4 , respectively; demonstrating this symmetry. The radial energy 
density of this solution is given in Figure 1 and one sees that it corresponds to a hollow 
ball. 



5 APPROXIMATE RADIALLY SYMMETRIC SKYRME FIELDS 



What about further genuine solutions? In general, our method does not give us further 
solutions but it is a matter of simple algebra to show that the condition ( |30D is true for 
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r 



Figure 1: Energy profile of the non-topological solution, 
any N > 2 when the modulus of the vector / is some power of (1 + |£| 2 ). In fact, we have 

N = 2, / = (!,£)*> (37) 

N = 3, f = (l,V2^ef, (38) 



N = n, / = (/o,/i,-..,/n-i) t : fi = ?yJC£?, (40) 

where denotes the binomial coefficients. Note that in this case, the constant K, in 

( p0|) is equal to the degree of the vector /: ie JC = n. 

Using the condition ([30]) the integrals involving P in the energy ( |I8|) can be evaluated 
analytically, 

J\f = n, 1={N-I) 2 = n 2 . (41) 

Using the analyticity of the projector P , it is straightforward to verify that the baryon 
number of this field is B = n, ie the degree of /. 



When ^ the energy (|18|) for the ansatz ( |15|) becomes 

I r gj 

E = — dr{A N g 2 r 2 + 2X sin 2 g(l + g 2 )+l - 

07T J 



sm 4 g 



2 2 

+m- r 



2<M , 1 ({N-l)2g 



(N - 1) (l - cos(jj-)) + 1 - cos 



N 



}. (42) 
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Minimising fl4*2| ) given (fll]) leads to the following equation for the profile function 

,2 



g„ [A N + 2n 



where A N is given by 



sin 2 g s 



+ 2A N — + — — ( n (p r - 1) - n — — 



— m 



9: 

/N — 1 
" V N 



r- 



r- 



2^ 

sin [ — ) + sin 



N 



0, (43) 



Solving fl43|) to determine g and then calculating the energy of the configuration we 
find that, for small m^, the energy for these configurations is a little higher than the 
energy of the SU(2) embedded ansatz with the same baryon number B when the mass is 
zero. However, when the mass increases, the picture changes. 



We have looked at field configurations corresponding to B = 2 — 4 for the SU(2) 
embeddings and for the SU(N) spherical symmetric fields (|38|)-(f4~0D where N = B + 1 (ie 
SU(3) for B = 2) and studied the dependence of their energies on m^. In all cases at low 
values of the mass the embeddings have lower energies while as the mass increases the 
energies increase. However, as the embedding energies increase faster for all low B there 
is a value of m n above which the embedding energy is higher. Unfortunately, this value 
of m n is quite large and it increases with the increase of B. 



These results are summarised in Table 1, which gives values of the energy for different 
values of the mass, and in Figure 2 where we present the dependence on m n of the energies 
for the embeddings and for the radially symmetric fields (|3~8|)-(fi(i|). Note that the energy 
per skyrmion of the harmonic ansatz configuration is always lower than the energy of a 
single skyrmion. 





SU(2) 


SU(2) 


SU(3) 


SU(2) 


SU(4) 


SU(2) 


SU(5) 


m n 


Eb=i 


Eb=2 


Eb=2 


Eb=3 


Eb=3 


Eb=3 


Eb=4 





1.232 


2.416 


2.444 


3.553 


3.644 


4.546 


4.838 


0.2 


1.247 


2.444 


2.472 


3.594 


3.683 


4.597 


4.886 


1 


1.416 


2.795 


2.808 


4.125 


4.172 


5.270 


5.520 


2.23 


1.693 


3.381 


3.370 


5.021 


5.006 


6.419 


6.615 


7 


2.510 


5.101 


5.030 


7.634 


7.478 


9.776 


9.880 


30 


4.783 


9.836 


9.633 


14.793 


14.339 


18.971 


18.948 



Table 1: Mass dependence of the energy for the radially symmetric configurations in the SU(2)- 
SU(5) models. 



Our configurations, like the exact solution of the SU(3) model mentioned above, all 
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Figure 2: Mass dependence of the SU(2) and SU(B + 1) harmonic map configurations for 
B = 2-3. 

have spherically symmetric energy density distributions (ie shell like structures). In Figure 
3 we present the curves of the energy density, as a function of the radius, for the field 
configurations mentioned above when m, = 0. We note that as the topological charge 
increases (and we consider the SU(N) model with larger N) the effective radius of the 
distribution also increases. 



6 SU(3) CASE 

In this section, we restrict our attention to the SU(3) model, take = and construct 
low energy states with baryon number from one up to six. From now on, N = 3 and so 
A N , given by (|T9|), becomes A N = 4/3. 



6.1 GENERAL DISCUSSION 

As in the previous section, we minimise (|I8D by first minimising the integrals M and I as 
functions of P and then minimising (|i~8|) with respect to the profile function g. 

Once again, Af is minimised by the so-called self-dual solutions of the Euclidean CP 2 
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r 



Figure 3: Energy profiles of the SU(B + 1) configurations for B = 2 — 6 for m n = 0. 

sigma model. They are given by (|25|) where / is any polynomial holomorphic vector /(£) 
and their energy Af is given by the degree of /. 



Next we note that the angular part of the baryon charge ([23|) coincides with the 
expression for the topological charge of the CP 2 sigma model and so simplifies to 



where dS = sin 9 d6 1 



2i(i + iei 2 )- 2 <*e<£ 



(44) 



To minimise for a configuration with a given baryon number 5, we take /(£) 
to be a holomorphic vector of degree S which, by construction, minimises N '. First we 
use the global SU(3) invariance of the model to reduce the number of parameters to the 
moduli space of the two dimensional sigma model, ie to 

( e + a B -i^ B - 1 + ■■■ + \ 



f 



(45) 



where all the coefficients are complex except which can be taken to be real. Then we 
substitute (0) for / of the form (f55| ) into I and minimise numerically the integral with 
respect to all the coefficients. Finally, treating Af = n and X as two fixed parameters, we 
minimize (0) by solving the resultant equation for g: 

N{gl - 1 



9r 



, . 3 sin q\ q r 3sin2o 
1 + 2AT +2 — + 



4 r z 



4 r , 



1 



sin 2 g 



0. 



(46) 



14 



An interesting feature of the SU(2) multi-skyrmion solutions is the shape of surfaces 
of constant energy or baryon density. In fact, the energy and the baryon densities of 
the skyrmion solutions look very similar. For the baryon density these surfaces look like 
hollow shell-like structures with holes in it, while for the energy densities the holes are 
partly filled in and so are represented by local minima ||. 

In order to investigate the situation for our SU(3) field configurations, we have to 
look at the components of / given in (^3) and study their effects on the density (pH[). 
Writing / = (K, L, M) 1 where K, L and M are polynomials of degree 5, B — 1 and B — 2 
respectively, the integrand of (|23| ) takes the form 



B = g r sin 2 s(l + |£| 2 ) 2 



\K e L - L^K\ 


2 + \- 






M^K\ 


2 + 


M S L - L ? M| 2 




1 2 + | 




2 + |M| 2 ) 2 



(47) 



Note that the integrand of (fT7|) is a scalar with respect to U (3) transformations applied 
to the vector /. Hence, any modifications of / which can be interpreted as such £7(3) 
transformations are symmetries of (|47|) . 

The radial factor g r sin 2 g in (H7I) indicates that if the angular part of the density 
vanishes, the baryon density will have radial holes going from the origin to infinity. For 
the density to vanish at some point we must require that the three factors in the numerator 
of (f47"D must vanish together, ie must have a common root. This is true, when the three 
polynomials R\ = K^L — L^K, R 2 = K^M — M^K and R3 = M^L — L^K have a common 
factor. However, these polynomials have 2(5 — 1), 25 — 3 and 2(5 — 2) roots, respectively; 
with, in addition, a possible root at infinity (ie the south pole of the sphere). By counting 
powers we see that the density does not vanish at £ = 00 unless L is a polynomial of 
degree less than 5 — 1. 

From this we conclude that the baryon density can have at most 25 — 3 holes but, in 
general, it is likely to have fewer holes if any. Of course, when some terms in (]47|) vanish, 
the expression may (but does not have to) have a local minimum. Note that this is in 
complete contrast with the SU(2) configurations of Houghton et al. |IJ which always have 
2(5 — 1) holes. In the 577(2) case, the vector / has only two components and so there is 
only one factor in the numerator of the baryon density which thus has 2(5 — 1) zeros. 
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6.2 SPECIFIC FIELDS 



In this section we present the detailed form of harmonic maps which are used in the 
construction of the 577(3) skyrmion field ansatze. 

First of all, the B = 1 discussed in section 5, is the SU(2) embedded skyrmion 

(ie the hedghog ansatz). Next we discuss field configurations for B = 2 — 6. In each case, 
having found the map which minimises X, we solve numerically (f46|) and determine the 
corresponding profile function g. In Figure 4 we present the energy profiles (as a function 
of r) of the resultant skyrmion field configurations. The profiles are given by the integrand 
of ( |I~8"D where the angular part of the energy, contained in Af and X, has been integrated. 
In Figure 5 we present the 9 angular dependence of the baryon densities for B = 2 — 4 
(no (j> dependence). 



30 
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20 
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5 



1 2 3 4 5 

r 

Figure 4: Radial energy profiles for B = 2 to B = 6. 

B = 2 

Using the ansatz (|4"5D, we have minimised X numerically and have found / to agree 
with the ansatz presented in section 5, ie to be given by 

/= (e,V2£, l)\ (48) 

For this field configuration |P + /| 2 /|/| 2 = 2/(1 + |£| 2 ) 2 and hence, as shown in Figure 5, 
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B = 2\ 



/ 



\ 
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\ 







Figure 5: The 8 angular dependence of the baryon density ([44|). 



the baryon and energy density are independent of the polar angles on the sphere. Thus 
the energy density of the B = 2 field represents a hollow sphere. 



The baryon density of this configuration is axially symmetric and has the shape of a torus 
with a sphere on top of it. In Figures 6a and 6b, we present plots of surfaces of constant 
baryon density for two different values. The values we have chosen are respectively 0.3 
and 0.7 times the maximum value of the topological density. (In all the graphs that follow, 
we always express the constant value for the curve as a fraction of the maximum density 
value). Notice that for low density value, the three skyrmion configuration has the shape 
of a pear, while for higher density values it looks like a ring under a small ball. 

The energy density has the same symmetry and has a virtually indistinguishable shape. 
This is also true for all the fields that we will present below. 

Note that as all components of / are monomials, a transformation £ — > £' = £e"*, for 
any a (ie a rotation around the 2- axis), can be interpreted as an SU(3) transformation. 
Hence the baryon density is invariant with respect to such transformations, ie it is axially 



B 



3 



The numerical minimisation of X leads to the following expression for /: 




(49) 
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3 <» 



Figure 6: Baryon density for B=3: a: level=0.3 b: level=0.7 

symmetric. 

Let us mention that the energy density of our configuration is remarkably similar to 
the density of a SU(2) configuration corresponding to three skyrmions in a mutually at- 
tractive channel [I2| and to the corresponding three monopole configuration [Tj|. Given 
the similarity of our SU(3) three skyrmion configuration to the equivalent SU(2) scatter- 
ing ones as well as to three monopoles one may expect that other monopole configurations 
which arise during the scattering process might also have their SU(3) analogues. Indeed, 
as we will see, this is the case for our SU(3) four skyrmion configuration. 

The baryon density for (^) does not vanish except when |£| 2 is infinite. This is the 
case as the three terms in the numerator of (|4"?| ) do not have common factors; however, 
as the second term of ([|Up is a polynomial of degree one, the baryon density vanishes for 
£ = oo. Indeed, we see in Figures 6a and 6b that the density vanishes on the negative 
part of the z-axis {9 = tt). 

B = A 
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Figure 7: Baryon density for B=4: a: level=0.4 b: level=0.7 



For B = 4 we find 



f=(e, 2.719U 2 , l)*. 



(50) 



This configuration also leads to the energy and baryon densities that are axially symmetric 
and they have the shape of two tori on top of each other. In Figures 7a and 7b, we present 
plots of the surfaces of the baryon densities at two values. 



Once again, the densities corresponding to fl5"0f) are invariant with respect to a rotation 
£ — > £' = £e* a . Note that the baryon density for (f49"D does vanish when £ is zero or when 
its modulus |£ 2 | is infinite. This happens since the three terms in the numerator of (|47|) 
have a single common factor at £ = 0, and the second term of / is a polynomial of degree 
two - implying once again that the baryon density vanishes when £ = oo. Indeed, this 
can be seen in Figures 7a and 7b; clearly the density vanishes along the z-axis (6 = and 
9 = 7i"). Once again, a similar configuration has been observed in the scattering of four 



SU (2) monopoles |R5| and skyrmions in an attractive channel |12 



B 
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The holomorphic vector for B = 5 is given by 

f = (e -2.7£,2£ 4 +l, 9/2£ 3 )'. (51) 

Note now that a transformation = (ie a 90° degree rotation around the z-axis) 

corresponds to a global SU(3) transformation. Hence the densities are invariant under 
such transformations. Let us add that the SU(2) embeddings have very different shapes 
and symmetries (in fact they are symmetric under 120° rotations). 

It is easy to check that the baryon density corresponding to the field in ( |5l| ) does not 
have any holes. Despite this, one can see holes in Figure 8a and 8b; they correspond to 
regions of low, but non-zero, baryon density values. 

Note that by taking / in the form close to (|51f) , ie 

/=(V + ^£, D? + C, Ee)\ (52) 
all the three terms in the numerator of fl4"?|) have zeros when 

e = % (53) 

which gives four holes in the baryon density. So, since our field (EH]) is not very different 
from ( p^D our densities have minima; corresponding to the holes ( |53"D partially filled in, 



by going from (p2| ) to (pl|). 
B = 6 



The holomorphic vector for B = 6 is given by 

/= (£ 6 + 3£, l-3£ 5 , ke)\ (54) 

where k was found, numerically, to be 7.06. Once again the baryon density of the field 
([54]) does not have any holes but has regions where it is small but non-zero (see Figures 
9a and 9b). These figures show that this configuration has an icosahedral symmetry and 
this leads us to the conclusion that, modulo an SU (3) global transformation, (|53J) must be 
invariant under the following transformation [14[]: £ — >• £' = £e* 27r//5 (ie a rotation by 72° 
around the z-axis); £ — > £' = — (which corresponds to 9 — » 7r — 9 and <p —>■ it — <fi) and 
£—>£' = (£ + b)(b£ + l)" 1 where b = 2cos(27r/5) = (v^— l)/2- This last transformation 
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Figure 8: Baryon density for B=5: a: level=0.4 b: level=0.6 




Figure 9: Baryon density for B=6: a: level=0.4 b: level=0.6 
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imposes a condition on k in ([54]): it is easy to see that the SU(3) transformation on / 
must be of the form U = R/ det(i?) 1 / 3 with 

/ 25 + 15a 10 + 5a (150 + lOO)^ 1 \ 



R 



10 + 5a 25 + 15a -(150 + 100 ajk' 1 
-£;(3 + 2a) fe(3 + 2a) 15 + 10a ; 



(55) 



where a = — (1 + y/E)/2. Imposing the condition that the rows and columns of R are 
orthogonal to each other implies that k = a/50 ~ 7.071 which is within the precision of 
our numerical minimisation programme. 

Having presented our field configurations we can now discuss some of their properties 
and compare them with the SU(2) embeddings. 

First of all, for B = 1 we have only the SU(2) embedding. Its energy and baryon 
density is in the shape of a ball. For B > 1 our field configurations are different from 
the SU(2) embeddings. However, the densities of the baryon densities for B < 4 are all 
axially symmetric (see Figure 5). The B = 2 configuration is radially symmetric and the 
baryon density corresponds to a shell (in contrast to the toroidal SU(2) one); the B = 3 
configuration corresponds to a single skyrmion located around the north pole of the S 2 
sphere and the other two are below the equator (spread out to form a charge two torus- 
like structure), while the B = 4 configuration consists of four baryons which are in the 
shape of two partially overlapping tori close to the equator of the sphere. The fields for 
B > 4 are more complicated, their baryon densities have fewer symmetries as seen from 
our figures. The baryon and energy densities for the case of B = 5 resemble a structure 
consisting of two deformed tori, close to the equator, with an additional ball at the north 
pole of the angular sphere while for B = 6 they form a structure which is icosahedrally 
symmetric. 

These shapes are very different from what was seen for SU(2) fields and, as we have 
discussed above, they have also different symmetries. 

In Figure 4, where we have plotted the energy profile functions for baryon numbers 
from two to six we note that the effective size of the baryons increases with the increasing 
baryon number - this is reflected in the shift to the right of the profile functions. 

In the following table we present the energy values of the resulting Skyrme fields. All 
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the numerical values of the energies are given in units of B and hence are close to unity. 
These values are then compared with the SU(2) skyrmion embeddings obtained using 
rational maps in [Q]. We see that both field configurations have similar values of energy, 
although the energies of the embeddings are marginally lower. 



B 


J(S77(3)) 


SU(3) En/Sk (Ansatz) 


SU(2) En/Sk (Ansatz) 


SU{2) En/Sk (Numerical) 


1 


1 


1.232 


1.232 


1.232 


2 


4 


1.222 


1.208 


1.171 


3 


10.65356 


1.215 


1.184 


1.143 


4 


18.04501 


1.184 


1.137 


1.116 


5 


27.26 


1.164 


1.147 


1.116 


6 


37.33 


1.1458 


1.137 


1.109 



Table 2: Energy of SU(3) harmonic ansatz compared to the energy of the SU{2) harmonic 
ansatz and the energy of the SU(2) solutions obtained numerically 0. 

7 CONCLUSIONS 

In this paper, we have discussed various static field configurations of the SU(N) Skyrme 
model. We have shown that, in addition to the obvious embeddings, any solution of 
the £77(2) model generates a solution of the S77(4) model. Unfortunately, this solution 
is topologically trivial (ie its baryon density vanishes identically) and its energy is four 
times the original SU(2) solution. 

Next we have generalised the harmonic map ansatz of Houghton at al |1J and showed 
that this ansatz has allowed us to find another exact solution, this time of, the SU(3) 
model. The baryon number of this solution is also zero and its energy density is radially 
symmetric. However, its total energy is less than four in topological units and we have 
argued that it represents a bound state of two skyrmions and two anti-skyrmions. 

Using our generalisation of the harmonic map ansatz we have then presented topolog- 
ically nontrivial field configurations of the SU(N) Skyrme model with radially symmetric 
energy densities. They correspond to B = N — 1 skyrmions in SU(N) models. In the 
massless case their energies have turned out also to be above those of the SU(2) embed- 
dings. However, when mass is added to the model, for sufficiently large masses, their 
energies can be be lower than the energies of the embeddings. 
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We have also looked at various field configurations of the SU (3) model. The energy 
and baryon densities of these SU(3) fields exhibit shell-like structures; in all cases, except 
for B = 1, they are different from the corresponding structures seen in the 577(2) model 
and more symmetrical. Their energies are slightly higher but comparable to those of the 
embeddings. Their different symmetry properties suggest to us that although these em- 
beddings have higher energies they may be reflections of real states of the model showing 
that the SU (3) model can have many interesting solutions. To see whether this expec- 
tation is correct one has to perform numerical simulations - this so far has not been 
done. 

Finally, our projector ansatz suggests that one might try to construct further ansatze 
involving two or more projectors. Such ansatze will then depend on more that one profile 
function. This topic is currently under investigation. 
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